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Abstract. Infrared (IR) luminosity of galaxies originating from dust thermal emission can be used as an indicator 
of the star formation rate (SFR). Inoue et al. II2000I IHK) have derived a formula for the conversion from dust IR 
luminosity to SFR by using the following three quantities: the fraction of Lyman continuum luminosity absorbed 
by gas (/), the fraction of UV luminosity absorbed by dust (e), and the fraction of dust heating from old (J^ 10* 
yr) stellar populations (77). We develop a method to estimate those three quantities based on the idea that the 
various way of SFR estimates from ultraviolet (UV) luminosity (2000 A luminosity). Ha luminosity, and dust IR 
luminosity should return the same SFR. After applying our method to samples of galaxies, the following results 
are obtained in our framework. First, our method is applied to a sample of star- forming galaxies, finding that 
/ ~ 0.6, e ~ 0.5, and 77 ~ 0.4 as representative values. Next, we apply the method to a starburst sample, which 
shows larger extinction than the star-forming galaxy sample. With the aid of /, e, and 77, we are able to estimate 
reliable SFRs from UV and/or IR luminosities. Moreover, the Hq luminosity, if the Ha extinction is corrected 
by using the Balmer decrement, is suitable for a statistical analysis of SFR, because the same correction factor 
for the Lyman continuum extinction (i.e., 1//) is applicable to both normal and starburst galaxies over all the 
range of SFR. The metallicity dependence of / and e is also tested: Only the latter proves to have a correlation 
with metallicity. As an extension of our result, the local {z = 0) comoving density of SFR can be estimated with 
our dust extinction corrections. We show that all UV, Ha, and IR comoving luminosity densities ai z = give a 
consistent SFR per comoving volume (~ 3 W~^h Mq yr~^ Mpc~^). Useful formulae for SFR estimate are listed. 

Key words. ISM: dust, extinction — galaxies: evolution — galaxies: ISM — galaxies: starburst — infrared: galaxies 
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1. Introduction 

During the history of the universe, galaxies have evolved 
forming stars. As a result, the present universe is filled 
with a large amount of stars and a significant amount of 
radiative energy originating from such stars. Therefore, 
tracing the star formation activity over all the history in 
the universe is fundamental to understand how the present 
universe has formed. To quantify the star formation activ- 
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ity. Star Formation Rate (SFR), defined as the stellar mass 
formed per unit time, is often estimated. 

The star formation activity, more specifically the SFR, 
can be traced by young (age ^ 10 ''-10* yr)^ stars. Since 
short-lived massive stars are certain to be produced most 
recently, the SFR is traced with the luminosity of massive 
stars. One of the well known and commonly used tracers 
of massive stars is Ha luminosity (e.g., Kennicutt 11983(1 . 
because Ha photons originate from the gas ionised by 
massive-star radiation. Since massive stars are the strong 
source for Ultra- Violet (UV) photons, UV luminosity is 
also used as an indicator of SFR. A general review for 
SFR indicators can be found in Kennicutt y^98a,). 



^ In this paper, the word "young" is used to specify the 
timescale on which the current SFR is traced. The word "old" 
is used otherwise. 
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In order to obtain a reliable estimate of the SFR, we 
have to take into account dust absorption (extinction)^ 
of UV and Ha light. The Ha luminosity (or more gen- 
erally, the luminosity of hydrogen recombination lines) is 
decreased also by the extinction of Lyman continuum pho- 
tons (e.g., Smith et al. lWTSl Inoue et aL lMiTl Inoue dOHTl 
Chariot & FaU'Min', Chariot & Longhetti EIMl Chariot 
et al. i2002. Dopita et al. 2003 ) . Quantifying the dust ab- 
sorption is generally crucial in deriving the star forma- 
tion history of a galaxy (e.g., Inoue et al. 120011 Hopkins 
et al. 120011 Kewley et al. 12002)1 or of the universe (e.g.. 
Floras et al lTM^ Steidel et al. ITU^ Meurer et al. 11^551) . 
Nevertheless, the correction for dust absorption usually re- 
quires an elaborate multi-wavelength or multi-band mod- 
elling (e.g., CalzettimnH). 

There are some SFR indicators that are virtually 
free from dust extinction. One of them is the infrared 
(IR) luminosity originating from dust continuum emission 
(Kennicutt '1998a). In this paper, we use the term IR to 
indicate the wavelength range where dust emission domi- 
nates the luminosity (^ 8-1000 /im). We call the total lu- 
minosity of dust emission "dust IR luminosity." Contrary 
to UV and Ha, dust IR luminosity traces the stellar radi- 
ation absorbed by dust. Thus, if a significant fraction of 
stellar radiation is absorbed by dust, detecting the dust IR 
emission is important. Recent ISO (e.g., Takeuchi et al. 
VWrn . COBE (e.g., Hauser & DweklWiTl . and SCUBA 
observations (e.g., Blain et al. '1999'; Barger et al. 2000) 
have shown a significant contribution of dust IR emission 
to the total light in the universe over a large part of the 
cosmic history. 

Theoretically the conversion factor between the dust 
IR luminosity and SFR is dependent on how efficiently 
stellar light is absorbed by dust and reprocessed in IR 
(e.g., Inoue et al. 120001 hereafter IHK). In the analytic 
formula of IHK, the conversion factor between IR lumi- 
nosity and SFR are described by using the following three 
parameters: /, e, and 77, where / is the fraction of Lyman 
continuum luminosity absorbed by hydrogen atoms, e is 
the nonionising photons from young stars absorbed by 
dust, and 77 is the fraction of IR luminosity originating 
from dust heating by old stars. 

Let us note some remarks about the three parame- 
ters. IHK assume that the ionizing photons do not escape 
out of galaxies; that is, the fraction 1 — / of the Lyman 
continuum luminosity is absorbed by dust grains. This 
can be justified for nearby galaxies, whose escape frac- 
tion of Lyman continuum photons is generally less than 
10% (e.g., Leitherer et al. [TM51 Deharveng et al. IM7T|l . 
For high-redshift (high-z) galaxies, it is still a matter of 
debate whether the escape fraction is high (Steidel et al. 
POOT) or low (Heckman et al. I^MTl Giallongo et al. IM?^ 
Fernandez-Soto et al. I2()()H(I . We can regard e as the frac- 
tion of UV (912-3650 A) photons absorbed by dust, partly 

^ Strictly speaking, the extinction is defined as the sum of 
scattering and absorption. In this paper, the term "dust ex- 
tinction" is used to indicate only the absorption. 



because almost all the radiative energy from young stars 
lies in the UV range, and partly because the dust absorp- 
tion is much more efficient in UV than in optical (Buat 
& Xu I1995|l . With respect to ry, we should keep in mind 
that the definition of 77 is temporal, not spatial as "cirrus". 
Indeed, young stars can heat dust outside Hii regions and 
it is shown that the UV heating is also efficient in all the 
interstellar medium (ISM) via a diffuse interstellar radia- 
tion field not spatially associated with Hll regions (Buat 
& XuT^, Waherbos & GreenawahTM^. Our definition 
of r] is also different from the cool dust fraction as defined 
in Lonsdale Persson & Helou H1987fl . because dust located 
out of Hii regions can be heated by UV radiation from 
massive stars and at the same time it can remain cool. 
Because of such a temporal definition of r/, we have to 
specify the timescale on which the current SFR is traced 
(Section EU. 

It is known that IR/UV flux ratio is a good indicator 
for the dust absorption in UV (Buat et al. 119991 Meurer 
et al. IMS* Witt & Cordon IMmH Panuzzo et al. ^M^ . 
Therefore, UV and IR luminosities are useful not only to 
derive SFR but also to correct UV luminosity for dust ab- 
sorption. By using Ha luminosity in addition to UV and 
IR luminosities, we develop a method to estimate /, e, 
and 77 (the principal quantities in the IHK formalism), and 
propose a way to obtain a reliable estimate of the SFR. 
Rosa-Gonzalez et al. H2002|l also treat e to estimate the 
SFR from dust IR luminosity, but we stress that / and 
rj are also important (IHK). Chariot et al. fi20lJ2| have 
compared the UV, IR, and Ha SFR estimators based on a 
fully consistent model, which includes different extinctions 
between young and old stellar populations, non-Balmer 
emission lines, and various types of star formation histo- 
ries (see also Chariot & FaU 20001) • Chariot et al. POU^ 
also use [O 11] line luminosity (see also Gallagher et al. 
^989), which is not included in our paper, to estimate 
SFR This paper adopts a different approach: instead of 
modelling the details of radiative processes, we develop 
an independent and simple way to extract the important 
quantities for SFR estimates, so that our model might be 
applied easily to large data sets. 

Hirashita et al. H2()()ll hereafter HOI) suggest that 
those three quantities, /, e, and rj, change as galaxies are 
enriched by metals. As a galaxy forms stars and recycles 
gas into ISM, the metallicity increases in some classes of 
models such as a closed-box model (e.g., Tinslev ri980|l . At 
the same time, dust grains are made from metals. In fact, 
metalhcity is related to dust content (e.g., Issa et al. 119901 
Schmidt & Boiler Lisenfeld & Ferrara lTMSI Dwek 

|B98, Hirashita 1999; Edmunds 2001). HOI consider that 
the metallicity evolution results in the increase of dust op- 
tical depth. Consequently / and e (and possibly ?/) can be 
affected by metallicity evolution. 

The aim of this paper is to develop a method to ob- 
servationally estimate the quantities important for deter- 
mining SFR. The luminosities in this paper are derived 
by assuming an isotropic radiation, which we expect to 
be reasonable in a statistical sense. First we reconstruct 



H. Hirashita et al.: Extinction and star formation rate 



the IHK formula to make it consistent with our treatment 
in this paper (Section [21). In Section |21 then, we explain 
how to derive the important quantities (/, e, and rj) for 
estimating SFR. In Sectional we present the samples to 
which we apply our method. The statistical properties of 
various SFR indicators are discussed in Section[Sl in which 
the metallicity dependence of / and e is also tested. In 
Section |B| we summarise our results and discuss possible 
application of our method to the cosmic Star Formation 
History (SFH) and to future survey data. 

2. Reconstruction of SFR formulae 

2.1. Reconsideration of the formula between Ljr and 
SFR 

We reconstruct the conversion formula from dust IR lu- 
minosity (Lir) to SFR following the method of IHK. [In 
Dale et al. H2()()l|l . this luminosity is called "Total IR lumi- 
nosity (TIR)."] A subset of the sample is commonly 
used to investigate the IR properties of galaxies, but the 
IRAS data are limited to the wavelength shorter than 120 
jjm. Therefore, the conversion formula from IRAS lumi- 
nosity defined between 40 and 120 /im to the total dust IR 
luminosity is useful. We call the luminosity in the 40-120 
jim range "FIR luminosity" (Lfir) , whose way of estimate 
can be seen in Lonsdale Persson & Helou ^WS7). Recently, 
data at longer wavelengths by ISO have enabled Dale et 
al. l!,2001 )) to estimate the ratio iiR/LpiR as a function 
of the IRAS 60 //m vs. 100 fim flux ratio. Hereafter we 
will use this their result unless otherwise stated. The ratio 
Lir/Z/fir is larger (> 2) for normal star-forming galaxies 
(such as spiral galaxies) than for starburst galaxies. While 
the past analyses based on starburst models (e.g., Meurer 
et al. I1999|l are not afi^ected by a small difference between 
LiR and Xfir (iiR/ipiR — 1-4; but see Calzetti et al. 
I2()()()l who derive im/ipiR — 1.8)'^, we have to be careful 
about the difference for normal galaxies. 

The current star formation activity of galaxies can be 
traced by quantifying the amount of young stars. Thus, 
the term "young" should indicate the timescale on which 
the current star formation is traced. This timescale is de- 
noted as tsF- For example, Inoue l|2()()2a|) adopts tsF ~ 10'' 
yr to trace the SFR. It is important to adopt a common 
timescale for all the SFR tracers. Since our aim is to work 
with monochromatic data near 2000 A, the best choice 
of the timescale is that appropriate for this wavelength. 
We choose isF = 10® yr unless otherwise stated because 
the 2000 A luminosity reaches its stationary value around 
10* yr in a constant SFR. Moreover, for large galaxies 
active in star formation but not necessarily starbursting, 
a constant SFR over 10® yr seems reasonable since the 
strong correlation found between the Ha and UV emis- 
sions (e.g., Buat et al. 120021 hereafter B02) argue for such 

^ For our lUE sample, which we treat as a starburst sam- 
ple, the application of Dale et al. 12UU1II to our data leads to 
Lm/Lpin — 2.0 on average. The uncertainty in Lir by 50% 
causes the change of /, e, and by ~ 20%. 
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a stationarity (but the stationarity for other samples is 
not necessarily supported; Sullivan et al. I2001|l . Such a 
hypothesis strongly simplifies the analysis and helps the 
understanding. Since the duration of the current star for- 
mation activity may be shorter in starburst galaxies, we 
also examine tsF — 10^ yr in Section 15.21 

IHK have established a procedure to derive the conver- 
sion formula which connects dust IR luminosity and SFR. 
They start from the following relation among luminosities 
originating from young stars (Petrosian et al. 11972( 1: 

^IR ~ I^Lya + (1 - /)iLyc + eL 

nonion ; 

(1) 

where L^^, ^Lya, -^Lyc, and Lnonion are various kinds of 
luminosities originating from young stars (luminosities of 
dust IR, Lya, Lyman continuum, and nonionising pho- 
tons, respectively), and / and e are the fraction of Lyman 
continuum luminosity absorbed by gas and the fraction of 
nonionising photon luminosity absorbed by dust, respec- 
tively. In this formalism, all the Lya photons are assumed 
to be absorbed by dust grains during some resonant scat- 
terings. This assumption can be justified even for dust- 
deficient galaxies with 1% of the Galactic dust-to-gas ra- 
tio (HOI), but we should be careful about this point if we 
apply our formula to objects from which a large amount 
of Lya photons leak for some reason. As we mentioned in 
Section ^ it is also assumed that the Lyman continuum 
photons are absorbed either by gas or by dust and do not 
escape out of the galaxy. 

All the three luminosities on the right-hand side in 
equation are related to the bolometric luminosity of 
the young stars, ^boi- order to obtain such relations, 
we have run Starburst99 (Leitherer et al. 11999^1 and made 
synthetic stellar spectra. We adopt the Salpeter stellar ini- 
tial mass function (IMF) with the upper and lower masses 
of 100 Mq and 0.1 Mq, respectively, constant SFR, and 
solar metallicity. We use the result at the age of 10® yr, 
and this timescale should be equal to tg-p. Then, we fi- 
nally obtain i^Lyc = O-lSL^^j and inonion = 0.87 L^l^ from 
the synthetic spectrum at tgp = 10® yr. In this paper, we 
adopt the same parameter set for the Starburst99 spec- 
trum unless otherwise stated. 

The Lya luminosity is estimated under Case B 
rOsterbrock iraS^ as (IHK) 

2 

Lhya — -^Ni^ycf'hl^hya , (2) 

where A^Lyc is the number of ionising photons emitted 
per unit time, hvi^ya is the energy of a Lya photon 
(1.63 X 10~^^ erg), and /' is the number fraction of Lyman 
continuum photons absorbed by gas (note that / is the lu- 
minosity fraction of Lyman continuum absorbed by gas). 
Although the relation between / and /' depends in a mi- 
nor way on the extinction law for Lyman continuum pho- 
tons and the shape of the Lyman continuum spectrum, 
/ ~ /' can be expected. There is a large difficulty in mod- 
elling the relation between / and /' because of the lack of 
knowledge about the extinction law for the Lyman contin- 
uum photons. Thus, we simply adopt / = /' throughout 
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this paper. By using the conversion from A^Lyc to -/^Lyc 
predicted by the synthesised spectrum of the Starburst99 
(^Lyc — 3.22 X lO^^iLyc in the cgs units), we obtain 
Lhya = 0.34/LLyc- Then, equation is reduced to 



rSF 



(0.13 - 0.085/ + 0.87e)Lg;;i. 



The Starburst99 resuh indicates that 
SFR 



Mq yr- 



1.79 X 10 



rSF 
-10 -'^bol 







(4) 



The dust IR luminosity, Lir, is the sum of two components 
originating from young (< tsp) stars and old (> tsp) stars 
(remember that tgp = 10* yr unless otherwise stated). 
Then the fraction of old stellar contribution to Ljr , rj, is 
defined as 



l!1 = {i-^)l 



IR, 



(5) 



Considering the energy balance of dust, rj is also inter- 
preted as the fraction of the energy input into dust (i.e., 
dust heating) from the old stellar population. Using equa- 
tions (O, (0)), and |(SJ), we obtain 



SFR _ 1.79 X lG-i"(l - ri) 
Mq yr-1 ~ 0.13 - 0.085/ -K0.87e 1^ 



(6) 



For the following convenience, we express the above con- 
version formula as 

SFR = aR(/, e, r;)LiR, (7) 
that is, the conversion factor Cir becomes 
1.79 X 10-i°(l -77) 



Cir(/, e, '7) 



[Mq yr-1 L-i] . (8) 



0.13 - 0.085/ -I- 0.87e 

We observe that the conversion factor is determined by a 
set of (/, e, 77). We should note that this numerical expres- 
sion for the conversion factor is based on the assumption 
that the star formation occurs at a constant rate for 10* 

yr- 

2.2. SFR from IR and UV emissions 

We formally define SFR(IR) as 

SFR(IR) =CitLiR, (9) 
where we define the conversion factor, Cf^, as follows: 
= Cir(/ = 0, e = 1, = 0) 

= 1.79 X 10^1" Mq yr-i Lq^ . (10) 

If the radiation field in a galaxy is dominated by the 
young stars and all the radiation from stars is once ab- 
sorbed by dust and recmitted in IR, SFR(IR) gives a 
reliable estimate for SFR (IHK). Indeed, Ljr = L^^j in 
such a case, and equation is the same as equation Q). 
This situation may be realised in starburst galaxies as as- 
sumed in Kennicutt t.l998b^ . Thus, we call the condition 
(/, e, 77) = (0, 1, 0) "dusty starburst approximation." Cf^ 
strongly depends on tgp (Appendix^). 



The intrinsic UV luminosity can be related to the SFR 
in a rather straightforward way. We adopt the 2000 A 
monochromatic luminosity ^2000 to trace UV (B02), and 
we express the SFR formula as 



(3) SFR(UV) = C2000 L 



2000 ■ 



(11) 



where C2000 can be calculated by using the Starburst99 
spectrum (with the same parameters as above in 
Section 12.11) without dust absorption: C2000 = 2.03 x 
10-40 (Mq yr-i)/(erg s'^ A'^). The advantage of us- 
ing 2000 A monochromatic luminosity is that (i) 2000 A 
is roughly the centre of the UV wavelength range (i.e., it 
traces the mean property of UV luminosity and extinc- 
tion) and (ii) a large number of UV data are available at 
2000 A. If there is no dust absorption, SFR(UV) gives a 
reliable estimate of SFR. 

Both SFR(IR) and SFR(UV) have their own disad- 
vantages. If we know 77, (1 — 77) SFR(IR) can be estimated 
from the observed Ljr. However, (1 — 77) SFR(IR) would 
systematically underestimate the SFR, because a part of 
the stellar radiation is not absorbed by dust. This under- 
estimate can be supplemented by SFR(UV) because most 
of the unabsorbed light from young stars is in UV. On the 
other hand, SFR(UV) can be formally estimated by mul- 
tiplying the observed 2000 A monochromatic luminosity 
with C200O1 but this always underestimates the SFR be- 
cause an appreciable amount of the radiation from young 
stars is absorbed by dust and reprocessed into IR. This 
underestimate can be supplemented by (1 — 77)SFR(IR). 
Therefore, the following sum of UV and IR SFRs is ex- 
pected to give a better approximation of the SFR: 



SFR ~ (1 - 77) SFR(IR) + SFR(UV) . 



(12) 



In the previous works such as Flores et al. I|1999|l and Buat 
et al. (Ungi, a simple sum SFR(IR) -I- SFR(UV) has been 
adopted. This type of simple sum is examined in Section 

The formula for SFR(IR) (eq.EJ is derived by assuming 
that all the stellar light is absorbed by dust. On the other 
hand, the conversion to SFR(UV) (eq. Ill|) is made by as- 
suming a theoretical stellar spectrum without extinction. 
Thus, we expect that equation (|12|l gives a reasonable es- 
timate for the SFR if there is no dust absorption or if the 
dust optical depth against the stellar light is significantly 
larger than 1. In the case where there is some dust absorp- 
tion, the shape of the absorbed spectrum in UV is modi- 
fied by the differential absorption, and the IR light comes 
from only a part of the UV radiative energy. Therefore, it 
is not obvious whether or not equation H12() is valid for an 
arbitrary value of dust extinction (e.g., E{B — U)). 

We check the validity of equation (|12|) . We start from 
the Starburst99 spectrum. Here, we only consider the con- 
tribution from young stars (i.e., 77 = 0), but we can apply 
the following consideration to 77 > by replacing SFR(IR) 
and LiR with (1 — 77) SFR(IR) and (1 — 77)Lir, respec- 
tively. Assuming the Calzetti extinction curve (Calzetti 
et al. I2000f) for the dust absorption, we extinguish the 
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Fig. 1. Star formation rates (SFRs) normalised with the 
true star formation rate. We show the part of the SFRs 
traced with UV and IR luminosities (dotted and dashed 
lines, respectively). The Calzetti extinction curve is as- 
sumed. The solid line shows the sum of the two SFRs and 
the correct SFR (normalised) should be 1. The small de- 
viation of the solid line from 1 means that the sum of UV 
and IR SFRs traces the true SFR. 

synthesised stellar spectrum as a function of E{B — V). 
The energy absorbed by dust is assumed to be equal 
to LiR. The energy absorbed by dust is estimated over 
all the range of the Starburst99 spectrum (100-1,600,000 
A)'* by using Calzetti et al.'s fitting formula. As a re- 
sult, we obtain £2000 (reduced according to the extinction) 
and LiR, and then SFR(IR) and SFR(UV) from equa- 
tions Q and (|ll|l . respectively. In this modelling, both 
SFR(IR) and SFR(UV) are proportional to SFR given 
in the Starburst99 calculation. Thus, wc should exam- 
ine if [SFR(IR) + SFR(UV)]/SFR = 1 is well satisfied. 
In Figure n we show SFR(IR)/SFR, SFR(UV)/SFR, and 
[SFR(IR) + SFR(UV)]/SFR as a function of E{B - V) 
(dashed, dotted, and sohd hues, respectively). The devi- 
ation of [SFR(IR) + SFR(UV)]/SFR from 1 is very small 
(^ 6%). 

We have used the Calzetti curve because it is the only 
extinction curve derived from a sample of galaxies. But 
its validity is checked only for starburst galaxies. However, 
even if we adopt the Milky Way extinction curve (Cardelli 
et al. 119891 with Ry ~ 3.1) instead of the Calzetti one, the 
difference of [SFR(IR) + SFR(UV)]/SFR from 1 is also 
^ 7%. Thus, we expect that the selection of a specific 
extinction curve does not change our conclusion that the 
sum of the IR and UV SFRs is an excellent indicator for 
the SFR. 

It may be worth mentioning qualitatively the follow- 
ing details, whose quantitative discussion nevertheless de- 
pends on the assumed stellar spectrum and extinction 

* Almost all (~ 90%) of the absorbed energy originates from 
the UV range. 



curve. The slight overestimate of SFR for E{B - U) ^ 0.1 
is due to a smaller extinction (i.e., a larger escape fraction) 
of the 2000 A light than that of the total UV hght. Because 
of this, the 2000 A monochromatic luminosity slightly 
overestimates the UV SFR. This point can be seen later in 
Figure 13 ( Section EHJ. On the contrary, SFR(IR) tends to 
underestimate the SFR even in the high-extinction limit 
as we can see for E(Ii — U) ^ 0.5. This is because a small 
part of energy escapes from galaxies in optical and near- 
infrared wavelengths. 

3. Estimates of the principal parameters 

We explain how to determine the parameters /, e, and 77, 
all of which appear in the IHK formula (eq. . We use the 
UV (2000 A), Ha, and dust IR luminosities to derive those 
three quantities. We illustrate our procedures in Figure 12 

3.1. e 

The parameter e is the fraction of UV photons absorbed 
by dust. Here, the UV wavelength range is defined between 
912 A and 3650 A. The intrinsic UV luminosity is denoted 
as Luv- The observational UV luminosity Luv can then 
be expressed as 

iuv = (1 - e) iuv ■ (13) 

Both Luv and Lyv defined as the luminosities in the 
UV wavelength range. 

As mentioned in Section |21 we use the 2000 A 
monochromatic luminosity to estimate the UV SFR. 
Therefore, we have to relate e to 62000 (defined as the 
fraction of 2000 A monochromatic luminosity absorbed 
by dust). This relation depends on assumed models. One 
of the simplest ways to avoid such model dependence is to 
assume £2000 — £• 

We quantify the uncertainty of this assumption £2000 = 
e. We adopt the Calzetti extinction curve to determine 
the absorbed fraction of UV light as a function of the 
extinction at 2000 A. For a certain value of the 2000 A 
extinction and the Starburst99 spectrum, we obtain the 
monochromatic luminosity after dust absorption for each 
wavelength. By integrating the monochromatic luminosity 
in all the UV wavelength range (912 - 3650 A) and com- 
paring this integration with the integration of the spec- 
trum before dust absorption, we obtain the fraction of 
UV light absorbed by dust. Consequently, we obtain e as 
a function of £2000- The e - £2000 relation is plotted in 
Figure 13 (sohd fine) . We find that £2000 tends to under- 
estimate £ but that the difference between the two is at 
most ~ 0.06. The difference in £ propagates to the esti- 
mates of / and T] in our framework described later, causing 
the uncertainty of at most ~ 20% in those two parame- 
ters. However, we should keep in mind that the conver- 
sion from £ to £2000 always has uncertainty if a specific 
model (i.e., a set of extinction curve and intrinsic spec- 
trum) is selected. For example, if the Galactic extinction 
curve (Cardelli et al. 1989; with Ry = 3.1) is assumed. 
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SFR(UV)+(l-r|)SFR(IR) = SFR(UV)/(l-8) = SFR(Ha)/f 

I I 
® © 



f, 8, r| Cm SFR(IR,IHK) 

Fig. 2. Summary of our method. We use three observational data (dust IR, UV, and Ha luminosities) and derive 
three parameters /, e, and 77, which are important for the conversion factor (Cir) from dust IR luminosity to star 
formation rate (SFR). The method is based on the consistency of the various SFRs derived from (1) the combination 
of UV and dust IR luminosities, (2) UV luminosity corrected for dust extinction, and (3) Ha luminosity corrected for 
dust extinction. Finally, using the three parameters, we derive Cir, which gives another independent estimate for SFR 
(denoted as SFR(IR, IHK)). 



we obtain the dashed line in Figure |31 In this case also, 
the difference from e — 62000 is small enough. We prefer to 
avoid the model dependence by simply adopting 62000 — £• 

The IR/UV flux ratio can be used to estimate 6 
(Meurer et al. \T^M Buat et al. \T^M B02). In principle, 
we can derive the extinction at 2000 A from a reasonable 
model of the radiative transfer such as Xu & Buat H 1995(1 if 
all the required data are available. Since the treatment of 
radiative transfer is beyond the scope of our simple anal- 
ysis, we adopt a simplified relation between IR/UV ratio 
and 2000 A extinction. For the SFG sample (see Section 0] 
for the sample description), we use the calibration of Buat 
et al. but expressed in F1R/F2000 (dust IR vs. 2000 

A) flux ratio instead of -FFiR/i^20oo (IRAS FIR vs. 2000 
A) one since the extrapolation of the total IR flux from 
that observed by IRAS has been made possible with the 
ISOPHOT observations (Dale et aL OUUTll . We adopt the 
same definition of i^A (A is a wavelength in UV in units of 
A) as that in Buat et al. H1999|l . i.e. Fx = A/a, where fx 



is the flux density per wavelength. Then we flnally derive 

A2000 - 0.21 + 0.75 log(FiR/F2ooo) 

+ 0.47[log(FiR/F2ooo)]' , (14) 

where A2000 (mag) is the extinction at 2000 A. We note 
that i^iR/F2ooo = iiR/(A2ooo ^2000), where A2000 = 2000 
A. This relation is derived from the energetics between IR 
and UV and is found rather insensitive to the details of the 
stellar spectrum if there is such an ongoing star formation 
activity as is seen in our samples (Buat & Xu 11996(1 . The 
fraction of 2000 A light absorbed by dust can be related 
to A2000 as 

6 - 62000 = 1 - lO-'^^o^/^.s _ (^5) 

We also investigate starburst galaxies observed by lUE 
as described in Section 0| For this lUE sample, we will 
use the extinction estimated by Calzetti et al. 1(2000(1 . The 
extinction at 1600 A can be observationally estimated for 
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0.8 
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0.4 

0.2 
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Galactic extinction 
e ~ £2000 



0.0 0.2 0.4 0.6 0.8 

£2000 



1.0 



Fig. 3. Relation between e (fraction of UV luminosity ab- 
sorbed by dust) and £2000 (fraction of 2000 A monochro- 
matic luminosity absorbed by dust) modelled for two types 
of extinction curves. The solid and dashed lines show the 
result for the Calzetti and Galactic extinction curves. The 
dotted line presents e = e2ooo- 



the lUE sample as (Calzetti EOHI 



A 



1600 



2.5 log 



1 Fu- 



1.68 F] 



1600 



(16) 



where i^ieoo is defined in the same way as ^20001 i-e., the 
flux density at 1600 A multiplied by 1600 A. For nearby 
galaxies, Figoo — -F2000 (Deharveng et al. 119941 Buat et 
al. I1999B . and thus we can use IR vs. 2000 A luminosity 
ratio to estimate ^iR/i^i6oo- Since A2000 = O.OAigoo under 
the Calzetti extinction curve, e2oooj which is assumed to 
be equal to e, is estimated by 



£2000 



= 1 - 10-"-'''4i6Oo/2-5 



(17) 



3.3. f 

We start from the following relation between SFR and 
Lyman continuum luminosity: 



SFR = Cl^^^ 



(20) 



where the Starburst99 spectrum indicates that CiLyc = 
3.45 X 10-"^ {Mq yr-i)/(erg s'^). We also obtain 



SFR = Cjvl.c ^'l 



(21) 



where A^Lyc is the number of the ionising photons emitted 
per unit time, and Cat^^^ = 1.07 x lO"^^ (M© yr-i)/s-i 
according to the Starburst99 spectrum. 

Since the Ha luminosity traces the amount of ionis- 
ing photons, we connect the above expressions to the Ha 
luminosity. According to Deharveng et al. H2001|l . 



/ATLy, = 7.34 X 10"L 



Ha 1 



(22) 



where is the Ha luminosity corrected for dust absorp- 
tion by using the Balmer decrement (i.e., observed Ha lu- 
minosity is multiplied by io^(Ha)/2.5^ where A(Ha) (mag) 
is the extinction for the Ha photons), and the quantities 
are expressed in the cgs units. This relation is determined 
from the Case B condition. Then we finally obtain the 
following relation from equations H21|l and (|22|l : 



SFR — Chq L'^a/ f ' 



(23) 



7.34 10" C 



7.89 10 ^'^ {Mq yr ^)/(erg s ^). The conversion factors 



where we obtain numerically Chq 

^0 yr" 

C^LycJ C'^VLycJ ^nd Chq are similar to the past estimates 
(e.g., Kennicutt ll998a|l . 

We estimate SFR in the left-hand side in equation H23|l 
from UV luminosity, dust IR luminosity, and -q (known by 
the method in Section f3.2|l by using equation (|12|l . By 
using equations and ((111) i^i addition, / is obtained: 



/ 



(1 



ri) Cf^ LiR -I- C2000 -^2000 



(24) 



3.2. T] 

We start from equation p2|) . where the basic idea is that 
the SFR measured with the IR emission is the SFR lost 
from the UV light because of the extinction. The left-hand 
side of equation l(T^ can be estimated by correcting the 
UV luminosity for dust absorption, i.e., 

SFR = SFR(UV)/(1 - e) . (18) 

For the right-hand side of equation (|12|l . SFR(IR) is es- 
timated from iiR (eq. El while SFR(UV) from L2000 (eq. 
Ill|) . Since e is known after applying the method in Section 
13.11 rj is determined from the following equation derived 
from equations iP, (HTJ, (II21 and lfTH|) : 

(1 - 77) Cf^ LiR = C2000 -^2000 ■ (19) 



4. Sample selection 

One of the samples of nearby star-forming galaxies whose 
extinction properties and star formation rates are well ex- 
amined is the SFG sample of B02. The sample consists 
of spiral and irregular galaxies located in clusters (Coma, 
Abell 1367, Cancer and Virgo). B02 also treat a starburst 
sample observed by lUE. We also use this sample {lUE 
sample) in order to test the applicability of our method 
and to investigate the difference in properties. We adopt 
only the galaxies with a good measurement of Balmer 
decrement (i.e., with a direct measurement of the under- 
lying stellar absorption). For the lUE sample, we only use 
the galaxies whose angular diameter is less than 1.5 ar- 
cmin in order to avoid the small aperture problem of lUE. 
The original quantities and the details are listed in Tables 
1 and 2 of B02 for the SFG sample and the lUE sample, 
respectively (see also the references therein). 
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The observational quantities used for our analysis are 
^Ha (t^^ Ha luminosity corrected for dust absorption by 
using the Balmer decrement), ^2000 (the monochromatic 
luminosity at 2000 A), and Lm (the dust IR luminosity). 
The dust IR luminosity is converted from the IRAS FIR 
luminosity by the method described in Dale et al. 1)2001(1 . 
The three quantities for the SFG sample and the lUE sam- 
ple are listed in TablesHandEJ respectively. We also exam- 
ine the dependence on metallicity, which is traced with the 
oxygen abundance. The oxygen abundance is taken from 
Gavazzi et al. (2003, in preparation) for the SFG sample 
and from Calzetti et al. H1994|l for the lUE sample. We 
call 12 -I- log(0/H) metallicity in this paper. The solar 
metallicity corresponds to 12 -I- log (0/H) — 8.93 (Anders 
& GrevesselMI Cox mi . 

From the arguments on dust temperature and equiv- 
alent widths of the Balmer lines, B02 show that the lUE 
sample consists of galaxies with higher star formation ac- 
tivity than the SFG sample. The lUE sample can be re- 
garded as a class of "starburst" galaxies, while the SFG 
sample can be representative of "normal" star-forming 
galaxies (spiral galaxies and irregular galaxies). We use 
the term "normal" and "starburst" to indicate a rough 
classification of star formation activity in this paper, and 
those two classes are often identified with the SFG sample 
and the lUE sample, respectively. 

5. Results 

In this section, we first apply our method to the SFG 
sample (Section I5.1|l . because the lUE sample may still 
have an effect of the small aperture even after we put 
the criterion for the angular size. The robustness of the 
method against the model is also examined for the SFG 
sample (Section I5.2|l . Then, we apply our method to the 
lUE sample in Section 15.31 The SFRs derived from the 
IHK method is compared with the best-estimate SFR in 
order to test the validity of the IHK formula fSection l5.4(l . 
Other SFR estimators are also examined in Section 15.51 
The metallicity dependence of some of the quantities (/ 
and e) is investigated in Section f5.6l in order to test the hy- 
pothesis in HOI. The luminosities and the SFR conversion 
factors are summarised in Tabled 

5.1. f, e, and rj for the SFG sample 

One of the characteristics of the SFG sample is that the 
Ha (corrected for dust absorption by the Balmer decre- 
ment) to UV flux ratio is lower than that expected for 
f — 1 under a constant SFR over 10^ yr (B02). This may 
indicate that some significant fraction of ionising photons 
is absorbed by dust grains. Inoue (j2001l . after analysing 
the individual H 11 regions of some Local Group galaxies, 
also reaches the same conclusion. Chariot et al. (|2002|l 
also concluded the same thing from their analysis of the 
Stromlo-APM redshift survey data. Thus, we expect to 
find a fraction / significantly lower than 1. 



In Tabled we list /, e, and 77 for each galaxy. By defini- 
tion, / should be between and 1, but only CGCG 119054 
shows / significantly larger than 1. We consider this to 
be due to the overcorrection of Ha absorption. Indeed, 
A(Ha) — 2.47 mag is the largest of all the SFG sample 
and the correction of dust absorption is as large as a factor 
of 9.7. With such a large extinction, the extinguished H/3 
fiux measurement could be very uncertain. As a result, 
the extinction derived from the Balmer decrement mea- 
surement could have a significant uncertainty. Therefore, 
we omit CGCG 119054 in the following analysis. 

The mean values and the standard deviations (ct) of /, 
e, and 77 are shown at the bottom of Tabled We calculated 
those values excluding CGCG 119054. As expected at the 
beginning of this subsection, / is significantly smaller than 
1. The mean value of / (0.57) indicates that about 40% of 
the ionising photons are directly absorbed by dust grains 
before being processed into recombination lines. This can 
be a reason for the systematic underestimate of Ha SFR 
relative to SFRs from other indicators (e.g.. Cram et al. 
1998: B02). We also observe from the mean of e (0.53) 
that the half of the UV is absorbed by dust and repro- 
cessed into IR. If we convert e = 0.53 ± 0.21 into A2000 
by using equation (|15|l and assuming e = £2000 (i-e., by 
A2000 = -2.51og(l - e)), we obtain A2000 = 0.82i°;^^. 
This confirms the result of B02. The value oi rj ^ 0.4 in- 
dicates that about 40% of the dust heating is due to old 
(age > tsF = 10* yr) stars, which have nothing to do with 
the current SFR. Misiriotis et al. (|2D0ip also find the con- 
tribution of old stellar populations to grain heating to be 
~ 40% for a sample of spiral galaxies. 

We also show the relations among the three quantities 
in Figure^! We see that there is no evidence for correlation 
either between / and e (r = 0.055; r is the correlation 
coefficient) or between / and rj {r — —0.017). There seems 
to be a tight relation between rj and e, but this tightness 
results from our formulation that determines both r/ and 
e as a fvmction of IR/UV flux ratio (see Sections 13. II and 
13. 2|) . In reality, there is a scatter in the relation between e 
and IR/UV flux ratio (see Fig. 1 of Buat et al. nM?jl . This 
scatter also disperse the rj-e relation and the standard 
deviation of rj would increase significantly. However, the 
mean values of 77 should still be ~ 0.4. 

5.2. Model uncertainty 

We examine how much the three quantities (/, e, and rj) 
are changed by the uncertainty in model assumptions. One 
of the sources for the uncertainty is the conversion between 
the /i?^S' FIR luminosity to the (total) dust IR luminosity. 
We have used Dale et al. ((2001(1 for this conversion. If we 
adopt Nagata et al. p002) for this conversion (their FIR2 
is used), the mean values of /, e, and 77 change to 0.62, 0.50, 
and 0.39, respectively (the standard deviations being 0.24, 
0.20, and 0.06, respectively). The difference in LiR/ipiR 
between Dale & Helou and Dale et al. ( 2002)) is 

smaller than that between Nagata et al. 1(20021:1 and Dale et 
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Name I/2000 L^o, Lm. 

erg erg erg s" 



VCC 25 


8, 


.91e+39 


2, 


.45e+41 


6, 


.76e+43 





.51 


0, 


,52 





.36 


VCC 66 


2, 


.34e+39 


1, 


.12e+41 


1, 


.95e+43 





.84 


0, 


,54 





.36 


VCC 89 


7, 


.94e+39 


2, 


.14e+41 


7, 


.08e+43 





.46 


0, 


,56 





.36 


VCC 92 


5, 


.37e+39 


2, 


.69e+41 


5, 


.Ole+43 





.83 


0, 


,57 





.37 


VCC 131 


1, 


.02e+39 


6, 


.46e+39 


4, 


.17e+42 





.16 


0, 


,36 





.39 


VCC 307 


8, 


.13e+39 


7, 


.24e+41 


2, 


.OOe+44 





.80 


0, 


,77 





.41 


VCC 318 


2, 


.04e+39 


3 


.39e+40 


4, 


.79e+42 





.50 


0, 


,22 





.48 


VCC 459 


6 


.31e+38 


1, 


.41e+40 


1, 


.17e+42 





.73 


0, 


,16 





.55 


VCC 664 


8, 


.51e+38 


1, 


.74e+40 


2, 


.24e+42 





.60 


0, 


,25 





.45 


VCC 692 


7, 


.94e+38 


1, 


.32e+40 


4, 


.37e+42 





.36 


0, 


,44 





.37 


VCC 801 


2, 


.4Qe+39 


1, 


.05e+41 


3 


.Q2e+43 





.61 


0, 


,64 





.38 


VCC 827 


1, 


.15e+39 


3, 


.31e+40 


2, 


.63e+43 





.27 


0, 


,76 





.40 


VCC 836 


1, 


.29e+39 


1, 


.20e+41 


3 


.89e+43 





.72 


0, 


,80 





.41 


VCC 938 


9 


.12e+38 


3 


.16e+40 


6, 


.03e+42 





.69 


0, 


,49 





.36 


VCC 1189 


6, 


.31e+38 


1, 


.OOe+40 


1, 


.62e+42 





.47 


0, 


,24 





.46 


VCC 1205 


1, 


.66e+39 


2, 


.19e+40 


8, 


.71e+42 





.29 


0, 


,43 





.37 


VCC 1379 


1, 


.82e+39 


2, 


.57e+40 


8, 


.13e+42 





.34 


0, 


,39 





.38 


VCC 1450 


1, 


.55e+39 


2, 


.40e+40 


7, 


.08e+42 





.36 


0, 


,39 





.38 


VCC 1554 


4, 


.07e+39 


2, 


.24e+41 


3 


.476+43 





.96 


0, 


,55 





.36 


VCC 1678 


6, 


.92e+38 


1, 


.29e+40 


1, 


.126+42 





.63 


0, 


,13 





.61 


CGCG 97087 


3 


.72e+40 


9, 


.33e+41 


2, 


.146+44 





.53 


0, 


,45 





.37 


CGCG 100004 


3, 


.80e+39 


1, 


.55e+41 


3 


.396+43 





.69 


0, 


,56 





.36 


CGCG 119029 


3, 


.16e+39 


2, 


.29e+41 


5, 


.756+43 





.80 


0, 


,71 





.39 


CGCG 119041 


5, 


.62e+38 


8, 


.13e+40 


5, 


.376+43 





.43 


0, 


,92 





.44 


CGCG 119043 


1, 


.48e+39 


1, 


.OOe+41 


3 


.636+43 





.61 


0, 


,77 





.41 


CGCG 119046 


4, 


.79e+39 


2, 


.19e+41 


2, 


.956+43 





.94 


0, 


,47 





.36 


CGCG 119047 


3 


.55e+39 


1, 


.95e+41 


7, 


.766+43 





.54 


0, 


,75 





.40 


CGCG 119054" 


1, 


.86e+39 


2, 


.34e+41 


2, 


.046+43 


1 


.90 


0, 


,61 





.37 


CGCG 119059 


1, 


.236+39 


4, 


.37e+40 


2, 


.756+43 





.34 


0, 


,75 





.40 


CGCG 160055 


1, 


.26e+40 


3, 


.39e+41 


1, 


.916+44 





.34 


0, 


,68 





.38 


CGCG 160067 


5, 


.62e+39 


3 


.09e+41 


5, 


.896+43 





.85 


0, 


,60 





.37 


CGCG 160139 


1, 


.02e+40 


2, 


.29e+41 


4, 


.176+43 





.55 


0, 


,36 





.39 


CGCG 160252 


5, 


.Ole+39 


3 


.80e+41 


1, 


.866+44 





.50 


0, 


,83 





.42 



Mean 0.57 0.53 0.40 

cr 0.21 0.21 0.06 



" This galaxy is not considered in taking the mean and a because / is larger than 1. 



al. 1)200 l|l . Therefore, with the available iiR/LpiR models, 
the parameters are determined within an uncertainty of 
- 10%. 

Rosa-Gonzalez et al. (|2002|l also calculated e. Some of 
their galaxies (13 galaxies) overlap with our lUE sample, 
whose e is derived later (Section 15.31 Table [21). They de- 
termined e in a similar manner as ours, but not the same 
(e.g., the way of the estimate of IR/UV flux ratio is differ- 
ent) . The agreement between our e and theirs is extremely 
good (the difference is within 0.05 except for Mrk 66, for 
which we obtain e — 0.69 while they derive e — 0.58). 
This supports the robustness of e for the lUE sample. 
For the SFG sample, the literature that analysed e is not 
found, and we cannot argue that other models support our 
derivation of e. 



We should remember that the conversion factors be- 
tween SFR and various luminosities depend on age, or 
more generally on SFH (e.g., Sullivan et al. I2001|l . In 
particular, the relative luminosity ratio between UV and 
Lyman continuum (or recombination lines) is sensitive to 
age. In a constant SFR, the luminosity of the Lyman con- 
tinuum reaches its stationary value at the age of ~ lO*" 
yr, while the UV luminosity becomes stationary at ^ 10* 
yr. In the above, we have adopted tgp — 10* yr for the 
age, but we could adopt tgp = 10'' yr to have an idea how 
much the quantities change in response to the SFH. Some 
of the conversion factors at 10'' yr change: C2000, C'lR, and 
Cf^ take the values described in Appendix^ while C^Lyc i 
CAfLyc, and Chq are the same as those at 10* yr. 

In the framework of Starburst99, it is very difficult to 
have a spectrum for an arbitrary SFH, but it is easy and 
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Table 2. lUE starburst sample. The sample is compiled in Table 2 of Buat et al. (2002). We only list and use galaxies 
whose angular size is less than 1.5 arcmin in order to avoid the small aperture effect. 



Name L2000 ^ho ^ir f ^ V 





erg s A 


erg s ^ 


erg s ^ 








Mrk 499 


8.51e+39 


4.57e+41 


2.00e+44 


0.32 


0.85 


-0.03 


Mrk 357 


6.46e+40 


2.09e+42 


3.89e+44 


0.50 


0.60 


-0.12 


IC 1586" 


4.37e+39 


7.76e+41 


8.71e+43 


1.20 


0.83 


-0.05 


Mrk 66 


4.90e+39 


1.02e+41 


4.37e+43 


0.25 


0.69 


-0.10 


NGC 5860" 


3.47e+39 


1.05e+42 


1.26e+44 


1.27 


0.89 


-0.00 


UGC 9560 


6.17e+38 


2.69e+40 


2.82e+42 


0.78 


0.54 


-0.13 


NGC 6090 


1.20e+40 


3.80e+42 


8.71e+44 


0.74 


0.94 


0.05 


IC 214 


1.02e+40 


1.15e+42 


1.05e+45 


0.20 


0.96 


0.08 


Tol 1924-416 


6.76e+39 


2.75e+41 


2.51e+43 


0.81 


0.49 


-0.14 


Haro 15 


1.38e+40 


3.02e+41 


1.23e+44 


0.26 


0.69 


-0.10 


NGC 6052 


5.50e+39 


3.89e+41 


2.75e+44 


0.23 


0.92 


0.02 


NGC 3125 


5.62e+38 


3.89e+40 


7.59e+42 


0.63 


0.77 


-0.07 


NGC 1510 


2.34e+38 


5.50e+39 


1.23e+42 


0.39 


0.57 


-0.12 


NGC 1614 


3.98e+39 


5.62e+42 


1.78e+45 


0.67 


0.99 


0.20 


NGC 7673 


5.01e+39 


4.37e+41 


1.17e+44 


0.52 


0.85 


-0.03 


NGC 7250 


8.71e+38 


2.88e+40 


1.07e+43 


0.32 


0.75 


-0.08 


NGC 5996 


1.38e+39 


1.74e+41 


5.37e+43 


0.50 


0.90 


0.00 


NGC 1140 


1.70e+39 


9.33e+40 


1.62e+43 


0.63 


0.70 


-0.09 


NGC 4194" 


2.40e+39 


1.91e+42 


3.39e+44 


1.04 


0.97 


0.11 


Mean 
a 








0.48 
0.20 


0.76 
0.15 


-0.04 
0.09 



" Those galaxies are not considered in taking the mean and a because / is larger than 1. 



Table 3. Luminosities and SFR conversion factors on the timescale of tsF = 10^ yr. 



Quantity 


Value 


Units 


Definition / comment 


1/2000 




erg s"^ A"^ 


2000 Amonochromatic 






erg s"^ 


corrected for Balmer decrement 


^IR 




org s^''" 


dust IR" 


C2OOO 


2.03 X 10"*" ^ 


(Mo yr-i)/(erg s'^ A'^) 


to be divided by (1 — e) 


Chq 


7.89 X 10""^ " 


(Me yr-i)/(erg s'^) 


to be divided by / 


^sb 


1.79 X 10"^" 


Me yr-^ 


dusty starburst approximation for IR 


Cm 


equation ||HJ 


Pie yr-' 


f,e, and 77 are necessary 



" The total luminosity of dust emission derived from the IRAS luminosity correcting for the contribution from longer (A > 
120 fim) wavelengths. Dale et al. Il2001t is used for the correction, but the difference between correction models is well less 
than 30% (e.g., the difference from Nagata et al. 12002^ . 
^ If tsF = 10'^ yr, the value becomes 3.18 x 10"*°. 
The value is unchanged for tsF ~ 10^ yr. 
If tsF = W'' yr, the value becomes 2.72 x 10"'°. 



computationally economical to change only the age under 
a constant SFH. Thus, we only change the age here. By us- 
ing another population synthesis code, we have also tested 
the difference between a constant SFR and an exponen- 
tially decaying SFR. As long as the exponential decaying 
timescale is comparable to (or longer than) the duration 
of the star formation (^sf), the change of th results is 
not so drastic as the difference between tsF = 10^ yr and 



tsF = 10 yr. If a galaxy has a decaying timescale much 
shorter than tsF, it would not be classified with a star- 
forming galaxy and would not included in our sample. 

If those coefficients for tsF — 10'' yr are adopted, we 
obtain the mean values / = 0.39, e = 0.53, and ij = 0.38 
(the standard deviations, cr, being 0.20, 0.21, and 0.06, 
respectively). It is natural that e does not change at all 
because it is determined from IR/UV ratio, which is inde- 
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Fig. 4. Diagrams of (a) / vs. e, (b) / vs. rj, and (c) rj vs. e for the SFG sample, where / is the fraction of ionising 
radiation absorbed by gas, e is the fraction of UV luminosity absorbed by dust, and 77 is the fraction of dust heating 
by the old (> 10^ yr) stellar population. 



pendent of the conversion factors. While 77 is not sensitive 
to the conversion factors, / changes significantly depend- 
ing on the assumed isF- Recalling that / is determined 
from the ratio between the SFR traced with Ha and the 
SFR traced with UV and IR, the sensitive change of / 
against tsF comes from the difference in the conversion 
factors for the UV and IR SFRs. Therefore, the typical 
age of the current star-forming activity is important if 
the age is shorter than 10® yr. We expect that our SFG 
sample has a continuous mode of the star formation be- 
cause the correlation between Ha and UV luminosities is 
good (B02). However, for a starburst sample such as the 
lUE sample, the duration of the present starburst could 
be shorter than 10® yr. Calzetti et al. ()1994|l assumed a 
constant SFR over 2 10^ yr for these galaxies. The corre- 
lation between Ha and UV luminosities is not so good as 
that for the SFG galaxies (B02), which also suggests that 
the sample has a diverse property in SFH on a timescale 
shorter than 10® yr. However, since we do not know the 
typical age of those sample, we adopt tsF = 10® yr also 



for this sample. If isF is shorter than 10® yr, / becomes 
smaller. Some mechanisms are proposed for a short-term 
variation of SFR on a galactic scale (Kamaya & Takeuchi 
I1H97I and references therein) . 

The spectral synthesis model used to derive the conver- 
sion factors is another source of uncertainty. We have seen 
above that the relative values between Cna and C2000 (or 
Cir) is the largest source for the uncertainty in /, although 
rj is determined quite robustly. Based on pegase synthesis 
code (Fioc & R.occa-Volmerange Iiy97|l . for example, the 
construction of the conversion factors are possible. Two of 
the conversion factors (Cua and C2000 in our notation) are 
thus obtained by Sullivan et al. 1)2001(1 for various age and 
metallicity. However, both of their C2000 and Chq are sim- 
ilar to ours (difference is < 20%). Therefore, the difference 
in / between the two synthesis codes is small and the dif- 
ference in tsF causes larger difference in /. The difference 
in metallicity affects Chq as largely as that in tsF , but the 
solar metallicity well approximates the metallicity of our 
samples. As long as the metallicity is between 1/10 and 3 
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solar metallicity, the major source for the uncertainty is 
the age, not the metalhcity. 

5.3. Application to the lUE sample 

The lUE sample is also analysed in the way described in 
Section ISTTl If we adopt the conversion factors for ^sf = 
10® yr, we obtain the results as listed in Table[21 For some 
of the galaxies, / is larger than 1, although / < 1 by 
definition. We suspect that this is due to the same reason 
as CGCG 119054 fSection l5.1|l . Indeed, the galaxies whose 
/ is larger than 1 have a large A(IIq;) (1.42, 1.69, and 1.96 
for IC 1586, NGC 5860, and NGC 4194, respectively). We 
exclude those galaxies from the following analysis. The 
mean values (±cr) are / = 0.48 (±0.20), e = 0.76 (±0.15), 
and rj = —0.04 (±0.09). If we apply the conversion factors 
at lO'' yr (i.e., we adopt tsF — 10^ yr), we obtain / — 0.30 
(±0.13), e = 0.76 (±0.15), and = -0.07 (±0.09). As 
seen in the SFG sample, / is sensitive to tsF also for the 
lUE sample. 

The extinctions for Lyman continuum and UV are 
both larger for the lUE sample than for the SFG sam- 
ple. For the lUE sample, about 50% of the Lyman con- 
tinuum photons and roughly 80% of the UV photons are 
absorbed by dust grains. Moreover almost all the dust 
heating source is the stellar population younger than 10® 
yr for the lUE sample, since 77 0. This is consistent with 
the starburst property of the lUE sample, whose current 
star formation activity dominates the luminosity of galax- 
ies. 

We examine the correlation between the quantities. 
There is no evidence for the correlations between / and e 
(r = —0.12) and between / and 77 (= —0.04). There seems 
to be a tight relation between e and 77, but this results 
from the same reason as the SFG sample fSection i5.1|l . In 
reality, the scatter should be much larger if we consider 
the scatter in the relation between IR/UV flux ratio and 
e. 

5.4. Validity of IHK 

One of our main aims is to examine the conversion formula 
of IHK (eq. [TJ. In order to see if the IHK formula works 
well or not, we should know the best estimate for the SFR 
first of all. We have shown that the SFR is traced very well 
by using IR and UV SFRs as equation p2l) . Therefore, 
in this paper the real SFR estimated observationally is 
defined as 

SFR(best) = (1 - 77) Qr, iiR + C^2ooo ^2000 • (25) 

Since /, 77 and e are supposed to be known at this 
step, we can examine the conversion factor for the IR 
SFR by using the IHK method. We define the following 
SFR(IR, IHK) by using the IHK conversion factor (eq.|Hl): 

SFR(IR, IHK) = Cir(/, e, 77) iiR . (26) 

Using the quantities (/, e, 77) derived for each galaxy in 
Tabled we show the relation between SFR(IR, IHK) and 



SFR(bcst) in FigurcEl We find that SFR(IR, IHK) agrees 
with SFR(best) within a difference of - 30% (for 70% of 
the SFG sample, the difference is within 10%). Therefore, 
if we know the three quantities, the IHK method approx- 
imates the SFR very well. 

For a general sample, we do not necessarily have all the 
three (UV, Ha, and dust IR) luminosities. In this case, 
Cir(/j v) cannot be obtained by our method. Thus, 
in the next subsection, we examine various SFRs derived 
from a limited number of luminosities. 

5.5. SFR 

Here we investigate SFRs derived from various indica- 
tors. In Figures Et-d, we compare SFR(UV), SFR(IR), 
SFR(IR,UV), and SFR(Ha) with SFR(best) for the SFG 
sample. Figure clearly shows that SFR(UV) under- 
estimates the SFR because of dust absorption. The ra- 
tio SFR(UV)/SFR(best) is equal to (1 - e) (eq. 
If we can estimate e (UV extinction) for each galaxy, 
SFR(UV)/(1 - e) gives the best estimate of SFR. 

The discrepancy between SFR(UV) and SFR(best) 
tends to be small for small SFR. For SFR(best) ^ 
1 Mq yr~^, SFR(UV) gives a good estimate for the real 
SFR (see also Bell & Kennicutt IMIT|l . This means that 
there is a positive correlation between SFR and e as 
pointed out also by Hopkins et al. H2001fl . This correlation 
may only reflect the size effect, since a large galaxy may 
tend to contain a lot of star-forming regions and at the 
same time a large optical depth of dust (Wang & Heckman 
[TMHI Buat & Burgarella lUMjl . 

We expect that SFR(UV)/(1 -e) with e = 0.53 gives a 
better estimate for the SFR. The dotted line in Figure 
EJi shows the relation SFR(UV)/(1 - e) = SFR(best) 
with e = 0.53. We observe that SFR(UV)/(1 - e) with 
e = 0.53 systematically overestimates the SFR for SFR ^ 
1 Mq yr~^, because the data points lie in the region 
SFR(UV)/(1 - e) > SFR(best). Thus, the dust correc- 
tion should be varied depending on the SFR. Because of 
a large variety in e, the scatter of SFR(UV) is larger than 
that of any other estimators. The dashed line in Figure 
El represents the relation applied to the lUE sample (i.e., 
e = 0.76). 

Figure EJs indicates that SFR(IR) estimates the SFR 
quite well for SFR ^ 1 Mq yr^^. However, we should 
keep in mind that this is the result of the cancellation of 
the following under- and overestimate (Kennicutt I1998al 
Inoue l2002a,l : SFR(IR) overestimates the SFR by a factor 
1/(1 — 77) because a part of the IR dust luminosity origi- 
nates from the old stellar population; SFR(IR) underesti- 
mates the SFR because a part of the radiation originating 
from young stars is not absorbed by dust, and thus is not 
traced by IR. However, SFR(IR) systematically underesti- 
mates SFR(best) when the SFR is lower than 1 Mq yr"-*^, 
because the major part of the energy is radiated in UV. 
Accordingly, SFR(UV) provides a reasonable estimate of 
SFR for SFR ^ 1 Mq yr"!. 
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Fig. 5. Star formation rate derived from the IHK method, SFR(IR, IHK), compared with the best-estimate SFR for 
the SFG sample and the lUE sample (a and b, respectively). 



The conversion factor for the IR luminosity can be 
tested by using the IHK conversion factor Cir(/, e, ry). In 
treating a sample of galaxies, a typical values of /, e, and 
T], for example the mean values (/ = 0.57, e = 0.53, and 
ry = 0.40 for the SFG sample), are useful. If we put those 
mean values, we find that Cir = 2.0 x 10~^° Mq yr^^ Lq^. 
This value is similar to CfJ^ = 1.79 x 10'^° Mq yr'^ L^^. 
If we use 2.0 x IQ-^" Mq yr^^ Lq^ instead of Cf^ to esti- 
mate SFR(IR), the data points in FigureEb shift upwards. 
In order not to complicate the figure, we shift the solid line 
down to the dotted line. The dashed line shows the same 
thing for Cir = 2.4 x 10^^° Mq yr~^ I^q^, which is repre- 
sentative for the lUE sample (Table0J. If /, e, and rj move 
their 1 a ranges, Cir = 1.3 - 3.5 x 10"^° Mq yr~ ^ Lq\ For 
a comparison, we should note that Buat & Xu p996|l de- 
rive a similar range Cir = 0.79 - 2.6 x 10"^° Mq yr'^ Lq\ 
where we assume -Lir/Lfir = 2.4 (the mean for the SFG 
sample) . 

We also examine the following SFR defined as a simple 
sum of IR and UV SFRs: 

SFR(IR, UV) = Cf^ LiR + C2000 ^2000 • (27) 

This kind of sum is adopted in Flores et al. H1999fl and 
Buat et al. (TSW). In Fi gure Et, we show the relation 
between SFR(IR, UV) and SFR(best). We observe that 
SFR(IR, UV) overestimates the SFR because the frac- 
tion 77 related to the old stars is not subtracted from 
LiR. However, the overestimate is not so large, ~ 60% 
at most. Moreover, the systematic SFR-dependent devia- 
tion, which is seen for SFR(UV) and SFR(IR), disappears 
by the combination of UV and IR SFRs. If we know a typ- 
ical 77 for a sample of galaxies, it is possible to statistically 
subtract the contribution from old stars by using equation 

The Ha SFR defined by the following expression is also 
tested: 

SFR(Ha) = Cho ^Ha ■ (28) 



SFR(Ha) is plotted against SFR(best) in Figure Hi. 
Because / is significantly smaller than 1, SFR(Ha) un- 
derestimates the SFR. The dotted line in Figure EJi shows 
the relation SFR(Ha)// = SFR(best) with / = 0.57 
(the mean value). This line reproduces the mean trend 
of the data over all the range of SFR. This trend strongly 
supports the usefulness of Ha luminosity as an indica- 
tor of SFR, because the Ha luminosity is independent 
of SFR(best), while SFR(best) includes dependence on 
SFR(UV) and SFR(IR) (thus the two SFRs in each of 
the figures a-c are not fully independent). Therefore, we 
conclude that SFR(Ha)// with / = 0.57 gives a good es- 
timate for SFR of the star-forming galaxies over the wide 
range of SFR. The dispersion in the figure is produced 
partly due to the different age in the present star for- 
mation activity, because Ha traces the star formation in 
recent 10^-10^ yr while UV and FIR traces all the star 
formation activity in recent ~ 10* yr or more. 

We should keep in mind that we have adopted the Ha 
luminosity corrected for j4(Ha). The analysis by Hopkins 
et al. H2001|l suggests that A(Ha) correlates with SFR 
(see also B02). Therefore, the conversion factor for ob- 
served Ha luminosity before the correction for A(Ha) is 
not universal but dependent on the SFR. It is also impor- 
tant that we should take into account not only ^(Ha) but 
also the Lyman continuum extinction in order to obtain a 
reliable estimate of the SFR. If we do not correct for the 
Lyman continuum extinction, the SFR is systematically 
underestimated by a factor of 2 (see also Inoue et al. 
IMIi . 

In Figure we examine the lUE sample. Since 77 < 
is not allowed by definition, we assume = if < 
0. This does not affect the following discussions because 
77 1 ^ 1 for the lUE sample. We use the conversion factor 
for tsF — 10* yr also for this sample. In Figure [7^, we 
show the hne SFR(UV)/(1 - e) = SFR(best) with e = 
0.53 (the mean value for the SFG sample) by the dotted 
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Fig. 6. Various SFRs of the SFG sample plotted against the best-estimate SFR. The solid lines show the equality 
of the two SFRs plotted in each figure, (a) SFR(UV) (without any extinction correction). The dotted and dashed 
lines indicate the "extinction-corrected" relations, SFR(UV)/(1 — e) — SFR(best) with e — 0.53 (the mean of the 
SFG sample) and with e = 0.76 (the mean of the lUE sample), respectively, (b) SFR(IR) converted directly from 
the dust IR luminosity by multiplying Cf^. The solid line shows SFR(IR) = SFR(best). The dotted and dashed lines 
present the trend for the SFG sample and the lUE sample, respectively, (c) SFR(IR, UV) = SFR(UV) -f SFR(IR). (d) 
SFR(Ha) derived from the Ha luminosity corrected for the Balmer decrement. The dotted and dashed lines show the 
relations corrected for the Lyman continuum extinction with / = 0.57 (the mean of the SFG sample) and / = 0.48 
(the mean of the lUE sample, respectively. 



line. Since the lUE sample is much obscured in UV, the 
correction with e — 0.53 underestimates the SFR. We also 
show SFR(UV)/(1 - e) = SFR(best) with e = 0.76 (the 
mean value for the lUE sample) by the dashed line, which 
fits the data points better. Therefore, when we correct the 
UV SFR of a galaxy for dust absorption, it is necessary 
to know if the galaxy is to be classified as a normal star- 
forming galaxy or a starburst. In this sense, there is no 
universal correction factor for the UV SFR. 

Figures[7|D shows that SFR(IR) approximates the SFR 
very well. This is because a large fraction of stellar light 
is absorbed by dust and reprocessed in IR in starburst 
galaxies. If we apply the mean values (/ = 0.48, e = 0.76, 
and 7? = 0), we obtain Cir = 2.4 x 10-^° Mq yr"! a 



value similar to Cf^ (i.e., SFR(IR, IHK)/SFR(IR) = 1.3). 
This is the reason why SFR(IR) gives a good estimate 
for the SFR of the starburst galaxies. The range expected 
from the 1 a variations of /, e, and 77 is Cir ~ 2.0-3.0 x 
10-10 Mq yr-1 L^K 

In Figured;, we plot SFR(IR,UV) against SFR(best) 
for the lUE sample. Since 77 ~ for the lUE sample, 
SFR(IR, UV) is almost equal to the best-estimate SFR 
(compare equations 1251 and I27() . Therefore, for starburst 
galaxies, the simple sum of the UV and IR SFRs gives the 
best estimate of SFR. 

The Ha SFR of the lUE sample is also shown in 
Figure di, in which the dashed line shows the relation 
SFR(Ha)// = SFR(best) with / = 0.48 (the mean value 
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Table 4. Recommended conversion factors from each luminosity to SFR on the timescale of isF = 10* yr. 



Luminosity^ 


Multiplyin 


g factor'" 


Comment 




normal'^ 


starburst'^ 




I'2000 


4.3 X 10~*° 


8.5 X lO"**" 


large dispersion in extinction 








systematic underestimate for ^ 1 Mq yr~' 




1.4 X 10"*^ 


1.6 X 10"^' 


similar factor for both ("universal") 








applicable to any SFR 


Lm. 


2.0 X 10"^° 


2.4 X 10"'" 


risk of underestimate for ^ 1 Mq yr~' 



Other formulae 

(1 - 77) C^R Lm + C2000 1/2000 r) necessary 

Cfjj LiR + C2000 L2000 systematic overestimate for normal galaxies 



^ is the Ha luminosity after the correction for the Ha extinction (yl(Ha)). Lm is the total luminosity of dust emission 

derived from IRAS 40-120 fim luminosity and 60 vs. 100 /im flux ratio (Dale et al. I200lt . 
^ The units are the same as in Table 

The SFG sample is assumed to be representative of normal star-forming galaxies. 
The lUE sample is assumed to be representative of starburst galaxies. 





Fig. 7. Same as Fig. |Hlbut for the lUE sample, li rj < 0, we adopt 7^ = 0. 
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for the lUE sample). The dotted hne presents the same re- 
lation with / — 0.57 (the mean for the SFG sample). The 
dashed line reproduces the mean trend of the data over 
all the range of SFR. Moreover, there is little difference 
between the dotted and dashed lines in Figure Et, which 
means that the same conversion factor can be used for Ha 
luminosity. Therefore, we suggest that SFR(Ha)// with 
/ ~ 0.5 gives a good estimate for SFR of both the SFG 
and lUE samples. Unlike the UV and IR conversion fac- 
tors, there is no systematic difference in the Ha conversion 
factor in all the SFR range. 

We summarise the above various ways of SFR estimate 
in Table 01 where we list the conversion factors applicable 
to both types of galaxies. We assume that the SFG sam- 
ple is representative of normal star-forming galaxies and 
that the lUE sample is representative of starburst galax- 
ies. The listed conversion factors are already corrected for 
dust effects (except for the Ha absorption, A(Ha), which 
should be estimated independently by Balmer decrement) 
and are our "recommended" values. 

5.6. Metallicity 

HOI suggest that / and e change as a function of metal- 
licity because the optical depth of dust can be related to 
dust-to-gas ratio (eqs. 11 and 13 of HOI) and dust-to-gas 
ratio increases as metallicity increases. We do not show the 
relation between rj and metallicity in this paper because 
it only shows rj ^ 0.4 for the SFG sample and rj ~ 0.0 for 
the lUE sample with a small scatter and without any cor- 
relation with metallicity. The scatter of rj should be larger 
in reality if we consider the scatter between the relation 
between IR/UV flux ratio and e (Section |OJ. 

In Figure |S1 we show the relation between (a) / and 
metallicity, and (b) e and metallicity for the SFG sam- 
ple. In Figure we do not find evidence for correlation 
(r = —0.13). There is a significant scatter over the range of 
< / < 1. However, Inoue et al. PUUT)) and Inoue (2001) 
show that there is a correlation between / and dust-to- 
gas ratio for Hii regions. They as well as HOI estimate 
the optical depth of dust for ionising photons by using the 
Stromgren sphere modelling of an H ii region. The result- 
ing optical depth becomes a function of the gas density 
of the Hii region and the ionising photon luminosity of 
the central star as well as dust-to-gas ratio CSpitzer ri978|l . 
Therefore, a possible interpretation on the large scatter of 
/ is that the gas density and/or the number of ionising 
photons per H ii region differ from galaxy to galaxy. HOI 
also assume that there is a tight relation between dust-to- 
gas ratio and metallicity. It is also shown that the scatter 
of dust-to-gas ratio can be so large that the correlation 
between dust-to-gas ratio and metallicity becomes weak 
(e.g., Hirashita et al. I2002|l . This would make the corre- 
lation between / and metallicity weak, even if there is a 
correlation between / and dust-to-gas ratio. Another rea- 
son for the large scatter is the presence of ionising photons 
escaping from H ii regions, which cannot be treated by the 



H II region modelling. The variety of geometry of dust dis- 
tribution (e.g., a dust-depleted region in the central parts 
of H II regions; Inoue l2002b|l can cause the large scatter of 
/• 

The correlation between e and metallicity (Figure 
r — 0.53) may directly support HOl's idea. HOI postulate 
a proportionality between the dust optical depth for the 
nonionising photons from young stars (Tnonion) and dust- 
to-gas mass ratio V as 

Tnonion = Tq (^J^^ , (29) 

where they calibrated the normalisation tq = 3.8 based 
on the Galactic condition. Since most of the nonionising 
photons are emitted in UV, Tnonion is related to e in a 
straightforward way: 

e = 1 - e~^"— . (30) 

By using the conversion from dust-to-gas ratio to metallic- 
ity as depicted in the solid line in HOl's Fig. 4, we finally 
obtain the model relation between metallicity and e. We 
show this relation in Figure|S|D by the solid line, which sig- 
nificantly overestimates the observed e. This implies that 
the normalisation tq is too large for the SFG sample. 

Then, we lower the normalisation tq to make the model 
applicable to the SFG sample. We adopt tq = 1.2 so that 
the mean e = 0.53 (i.e., Tnonion = 0.76) is satisfied at the 
mean metallicity (12-1- log (0/H) = 8.8). The e-metallicity 
relation under this lower normalisation is shown by the 
dotted line. This reproduce the observed trend quite well. 
This implies that HOl's model with tq = 1.2 can be appli- 
cable to star-forming galaxies. 

We examine the same relations for the lUE sample in 
Figure El Also for this sample, there is not any clear trend 
in the /-metallicity diagram (r — —0.07), but there is 
a correlation between e and metallicity (r — 0.74). Thus, 
we have confirmed the correlation between metallicity and 
UV extinction for the lUE sample (Heckman et al. ll998|l . 
The solid line in Figure Eb shows the prediction by HOI 
(i.e.. To = 3.8). Contrary to the SFG sample, tq = 3.8 is 
too small for the lUE sample. If we assume tq = 5.5 to 
satisfy the mean e — 0.78 (i.e., Tnonion = 1-3) at the mean 
metallicity (12 + log (0/H) = 8.6), we obtain the dashed 
line. However, even in this case, the data points cannot be 
reproduced because the extinction is extremely large even 
for low-metallicity galaxies in this sample. 

Therefore, although HOl's idea that there should be a 
relation between metallicity and extinction could be partly 
supported, the extinction is not described solely by a func- 
tion of metallicity as their original idea. The extinction is 
largely dependent on whether a galaxy is a "starburst" 
galaxy or a mild star-forming galaxy. The larger UV opti- 
cal depth for the lUE starburst galaxies implies that the 
star-forming regions of starburst galaxies are deeply em- 
bedded in dusty gas. Those two classes of galaxies may 
also be difi'erent in the IR/UV vs. UV spectral slope re- 
lation (e.g.. Bell 1200^ . which also implies a fundamental 
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Fig. 8. (a) / vs. metallicity, and (b) e vs. metallicity for the SFG sample. The sohd hne in (b) shows the model by 
HOI. The dotted line represents the result with our new normalisation for the optical depth of the SFG sample. 




Fig. 9. Same as Fig. |Hlbut for the lUE sample. The solid line shows the model by HOI. The dotted line represents 
the result with our new normalisation for the optical depth of the. 



difference in the extinction properties. Since our method 
is aimed at a simple treatment to allow for easy applica- 
tions, the physical modelling of the difference is beyond 
the scope of this paper and is left for future works. 

Moreover, there is an appreciable scatter in the e- 
metallicity relation. The scatter can be caused also by the 
varieties of following quantities: the inclination, the ge- 
ometry of dust distribution, the dispersion in the relation 
between dust-to-gas ratio and metallicity, etc. 

It is important that within each type of galaxies, there 
is a correlation between extinction and metallicity. This 
correlation is equivalent to the correlation between IR/UV 
flux ratio and metallicity. It is well known that there is 
a correlation between galaxy mass (or luminosity) and 
metallicity (Zaritsky et al. [TUM Richer & McCalli995, 
Garnett et al. I1997|l . There is also a correlation between 
mass (or luminosity) and IR/UV flux ratio (or extinction) 



(Wang & Heckman Heckman et al. [TM5I Buat & 

Burgarella UnnHI Buat et al.^M B02). Those correla- 
tions suggest that large galaxies work as larger reservoirs 
of gas and metals (and dust) (Wang & Heckman 11996(1 . 

6. Summary and discussion 

6.1. Summary 

In this paper, we analysed various SFR indicators (UV, 
IR, and Ha luminosities). Especially, we focused on the 
IHK formula that converts dust IR luminosity into SFR. 
For this conversion, the following three quantities are cru- 
cial: the fraction of ionising radiation absorbed by gas (/), 
the fraction of UV luminosity absorbed by dust (e), and 
the fraction of "old" stellar contribution to the total dust 
IR luminosity (77). Those three quantities were observa- 
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tionally estimated from the 2000 A monochromatic lumi- 
nosity, the Ha luminosity, and the dust IR luminosity for 
the SFG sample and the lUE sample compiled in B02. 

The SFG sample proved to have / = 0.57 ± 0.21, e ^ 
0.53±0.21, and 77 = 0.40±0.06. Those values mean that (1) 
about 40% of the ionising photons are directly absorbed 
by dust; (2) roughly half of the UV photons are absorbed 
by dust; (3) about 40% of the heating of dust is due to 
stars older than 10* yr. For the lUE sample, we found 
that / = 0.48 ± 0.20, e = 0.76 ± 0.15, and r] = -0.04 ± 
0.09. Therefore, the typical properties of the lUE sample 
is as follows: (1) about 50% of the ionising photons are 
absorbed by dust; (2) most ('^ 80%) of the UV photons 
are absorbed by dust; (3) almost all the heating source for 
dust grains is the stars younger than 10* yr. 

Based on those parameters, we examined the IHK for- 
mula. The SFR derived from this formula agrees almost 
exactly with the best-estimate SFR given by the combi- 
nation of IR and UV luminosities (eq. I25|l . This demon- 
strates the reliability of IHK's formula over a wide range 
in SFR (- 0.1-100 Mq yrj^ IHK's formula is different 
from that of Kcnnicutt l|1998b|l . where it is assumed that 
the dust IR luminosity is equal to the bolometric luminos- 
ity of young stars. This assumption is equivalent to the 
case of / = 0, e = 1, and 77 = 0. We call this assumption " 
dusty starburst approximation" . For the dusty starburst 
approximation, our Starburst99 calculation indicates the 
conversion factor of Cf^ = 1.79 x 10"!° Mq yr"! Lq^ 
Our result for the SFG sample implies that / ~ 0.57, 
e ~ 0.53, and 77 ~ 0.40 are applicable for nearby normal 
star-forming galaxies as a first approximation, and we ob- 
tain the conversion factor Cir ~ 2.0x lO^^*' Mq yr~^ ^o^y 
a value similar to that under the dusty starburst ap- 
proximation. This similarity comes from the two offset- 
ting effects as stated in Kennicutt (.1998a)): the contri- 
bution of old stars to the total IR luminosity and the 
escape of UV photons without being absorbed by dust 
(Section I5.5|l . The IHK formula works also for the lUE 
starburst sample with / ~ 0.48, e ~ 0.76, and ry ~ 0.0 
(Cir = 2.4 X 10-10 Mq yr'i Lq^). 

The SFG sample can be regarded as "normal" star- 
forming galaxies, and the lUE sample can be represen- 
tative of the "starburst" galaxies. After analysing various 
SFR indicators, we found the following (see also Table 0)): 

- UV: The UV SFR should be corrected for dust ex- 
tinction by multiplying 1/(1 — e). The correction factor 
depends largely on the property of individual galaxies, 
especially on the starburst/normal category. Among 
each population (especially among normal galaxies), e 
is systematically small for SFR ^ 1 Mq yr~^. Thus, 
extinction estimate for each galaxy by using e.g., IR 
luminosity, is important to obtain a reliable SFR. 
Panuzzo et al. H2003|l conclude that the UV luminos- 
ity corrected by using the IR/UV ratio is a reliable 
indicator of SFR. 

— Ha: If the Balmer decrement is measured precisely 
enough to correct for the extinction of Ha photons. 



Ha luminosity is the most "secure" estimator of SFR, 
This is partly because the correction factor (1//) for 
the Lyman continuum photons does not differ be- 
tween normal and starburst galaxies, and partly be- 
cause the there is no systematic trend of / with re- 
spect to the SFR. The dispersion of SFR(Ha) relative 
to the SFR(best) (SFR estimated from UV and IR) 
can be produced by age variation of the present star 
formation activity. 

— dust IR: The IR luminosity traces the SFR quite well. 
The conversion factor derived under the dusty star- 
burst approximation is applicable to both normal and 
starburst galaxies. There is a risk that the SFR is un- 
derestimated for SFR < 1 Mq yr-^ The IHK formula 
also provides us with a way to estimate the conversion 
factor if we know typical values of /, e, and 77. 

— Combination of IR and UV: The simple sum of IR 
and UV SFRs systematically overestimates the SFR 
of normal galaxies, because some fraction (7;) of IR lu- 
minosity is not related to recent star formation. If we 
know the typical 77, we can use equation 125|l to sub- 
tract the contribution from old stars to IR luminosity. 
This SFR is free from any SFR-dependent systemat- 
ics. If we know 77 for each galaxy, we obtain a reliable 
estimate of the SFR. 

The metallicity dependence of / and e was also tested. 
We found a correlation between e and metallicity for 
both samples, but we did not find any trend of the /- 
metallicity relation. The e-metallicity relation of the SFG 
sample implies lower extinction than that suggested by 
HOI (rnonion — 0.76 ou avcragc) . On the contrary, the 
lUE sample showed a higher extinction by 2.1 times 
(inonion — 1-6 On average). Compared at the same metal- 
licity level, the lUE sample has the UV optical depth 4.6 
times larger than the SFG sample. This is consistent with 
the picture that starburst galaxies are highly obscured by 
dust grains (e.g., Heckman et al. 11998(1 . 

6.2. Application to the cosmic SFH 

We comment on the application of our method to a cos- 
mological context. The cosmological evolution of galaxies 
is one of the main topics in the cosmic structure formation 
(e.g.. White & R.ees 11971^ . In particular, it has been an 
important and unsolved question how and when galaxies 
have formed stars (Tinsley & Danly 1980). Such a cosmic 
SFH as observationally derived by Madau et al. H1996|l 
provides some keys for the statistical view of galaxy evo- 
lution. 

Takeuchi et al. H2001|l apply IHK's formula to derive 
the cosmic SFH from IR data. They mainly use the num- 
ber counts by ISO and the cosmic IR background by 
COBE to constrain the comoving IR luminosity evolution. 
They multiply Cir (conversion factor of IHK) to the co- 
moving IR luminosity density and derive the cosmic SFH. 
They also show that if there is difference in / and e be- 
tween nearby and distant galaxies, the SFH derived for 
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the high-z {z ^ 1) universe is uncertain by a factor of 
~ 4. Therefore, the apphcation of our method to obtain a 
typical values for /, e, and rj for high-z galaxies is an in- 
teresting future topic. High-z galaxies might show a large 
dust extinction (i.e., small / and large e) (e.g., Heckman 
et al. Meurer et al. ITM^ Massarotti et al. I^IMIl or 

perhaps a small extinction (e.g., HOI). 

The typical properties of a local sample can be applied 
to infer the comoving density of SFR at z = 0. The lumi- 
nosity of galaxies per unit comoving volume at z = has 
been derived in a lot of literatures. In particular, the 2000 
A monochromatic luminosity and the IR luminosity can 
be converted to SFR by using the formula described in 
this paper. For 2000 A monochromatic luminosity, Buat 
et al. derive the following value for the comoving 

density at z = from the measurement of Treyer et al. 
ifT^ at z ~ 0.2: 

P20oo(2 = 0) = 8.9 ± 3.9 erg s^^A"^ Mpc"^ , (31) 

where h is the Hubble constant at z = normalised by 
100 km s~^ Mpc^^. This can be converted to the comoving 
SFR by multiplying 1/(1 — e) for extinction correction and 
C2000 for conversion if we know the luminosity-weighted 
mean of 1/(1 — e) for all the nearby galaxies. We can ten- 
tatively apply the mean e derived for the SFG sample 
(~ 0.5), because the nearby UV extinction is suggested to 
be smaller than the lUE sample and more similar to that 
of the SFG sample (Buat et al. [TM??| . Adopting e = 0.5 for 
the extinction correction, we obtain the comoving density 
of SFR at z = as 

Psfr(^ 0) = 3.6 ± 1.6 IQ-^h A/© yr^^ Mpc"^ . (32) 

In fact, we need a complete sample observed in both UV 
and IR to estimate dust extinction. Because the available 
UV samples are not large enough, we have to wait for 
future UV observations such as GALEX^ . 

The above e(= 0.5) might be overestimated for a 
UV selected sample (K. Xu, private communication). 
Contrary to it, an analysis of a currently available UV- 
selected sample by Sullivan et al. (Ij2001j^ proves a mean 
UV extinction to be 1.3 mag, larger than the value which 
we have adopted above (0.82 mag). Nevertheless their cal- 
culations are made using the Balmer decrement and the 
Calzetti extinction curve, and if their galaxies are similar 
to the SFG sample they probably overestimate the extinc- 
tion (e.g., B02). Thus, future observations are crucial to 
correct the comoving UV SFR for dust extinction even in 
the local universe. 

We can discuss the comoving SFR from IR data. 
Saunders et al. H1990|l estimate the comoving density of 
FIR at z = (see also Takeuchi et al. lM75)l : 

P¥m{z = 0) = 5.6 ± 0.6 lO^/i Mpc"^ . (33) 

If we multiply this with 2.4 (the mean value for the SFG 
sample) to obtain the total dust IR luminosity, we obtain 

^ http:/ /www. srl.caltech.edu/galex/ 



the comoving dust IR luminosity as 

Pik{z = 0) = 1.3 ± 0.1 lO^/i Lq Mpc-3 . (34) 

If we adopt Cir(/ = 0.57, e = 0.53, = 0.40) = 2.0 x 
10"^" (-^0 yr~^ ^0^) for the conversion factor from IR 
luminosity to SFR (typical for the SFG sample) , we obtain 
the following local comoving SFR density: 

Psfk{z = 0) = 2.7 ± 0.3 IQ-'^h M© yr^^ Mpc'^ , (35) 

in good agreement with equation (|32|l . An advantage of 
using IR luminosity is that we can apply a similar con- 
version factor whether a galaxy might be a normal star- 
forming one or a starburst one (Table However, we 
have to be careful about the systematic underestimate for 
SFR ^ 1 Mq yr""'^. If such low-SFR galaxies dominate 
the star formation activity in the local universe, the above 
SFR density is an underestimate. 

We can make the same kind of argument for the Ha 
comoving density derived by Gallego et al. I|1995|l (see also 
Tresse et al.EO^: 

/5Ha(2 = 0) = 2.5tj j 10^9/1 erg s'^ Mpc'^ . (36) 

In order to convert this to the comoving SFR density, we 
have to multiply Cho// (see eauation l23ll . Then we obtain 
PSFR = 2.0lJ:^ 10-V(/) h Mq yr-1 Mpc^^ where (/) is 
the typical / we should apply to the sample in Gallego et 
al. H1995|l . Since we do not know (/), we assume that it is 
equal to the mean / (0.57) in the SFG sample. Then, we 
obtain the following comoving SFR density: 

PSFRiz = 0) = S.btli 10"'/i Mq yr-i Mpc-^ . (37) 

Again we obtain a similar value as the above two estimates 
(equations 021 and 13 5|) . Those comoving SFRs also agree 
with Buat et al. |1999) . 

We also try to estimate the SFR by using both IR and 
UV data. According to equation (|25|) . the comoving SFR 
can be estimated as 

PsFuiz = 0) = (1 - ry)Cf^piR(z = 0) + C2ooo/O2ooo(^ = 0) 
= 3.2 ± 0.9 IQ-^h Mq yr-i Mpc^^ , (38) 

where we assume 77 — 0.4 (the mean value for the SFG 
sample). If we apply the simple sum of IR and UV SFRs 
(i.e., 77 = in the above estimate), we obtain a larger 
comoving SFR: psFR = (4.1±1.0) lO'^/i M© yr"! Mpc'^. 
Such a simple sum has been adopted by some authors 
(e.g., Buat et al. HHSi Flores et alT^UUI, but this may 
be an overestimate if a significant fraction of the dust IR 
luminosity originates from an old stellar population. 

6.3. Prospects for IR and UV surveys 

In the near future, a large IR sample will be obtained 
by SIRTF^, ASTRO-F^ and SOFIA « with the typical 

® http://sirtf.caltech.edu/ 

http:/ /www. ir.isas. ac.jp/ ASTRO-F/index-e. html 
* http://sofia.arc.nasa.gov/ 
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rcdshift z < 1. Even at the first step of the data release, 
the galaxy number count can be estimated, to which the 
modelling by e.g., Takeuchi et al. 1)2001(1 can be applied so 
as to obtain the comoving density of dust IR luminosity 
as a function of z (see also e.g., Gispert et al. ,2000j . If 
we only have dust IR luminosity, we can assume a typical 
values to apply the IHK formula (/ ~ 0.6, e ^ 0.5, and 
77 ~ 0.4) as a first approximation. 

More sensitive and high-resolution IR (or sub- mm) ob- 
servations by Herschel^, ALMA^o, SPICA^^ , etc. will de- 
tect a large number of high- 2: galaxies. For high- 2; galax- 
ies, the typical values for /, e, and i] derived in this paper 
may not be applicable because of the difference in dust 
amount, age, typical size, etc. HOI consider that the low- 
metallicity condition at high z makes the dust extinction 
less efficient than at low z. Then they suggest a high Cir 
at high z. Takeuchi et al. I|2001i,) also propose that if Cir is 
systematically larger at high z, we obtain a flat (or nearly 
constant) SFH from z ~ 1 to 2: ~ 5. On the contrary, some 
other works have suggested the importance of dust extinc- 
tion for UV photons for high-z galaxies (e.g., Meurer et 
al. n^^ Steidel et al. flM^ . More recently, Papovich et 
al. 1)2001(1 and Seibert et al. 1)2002(1 have shown that the 
correction factor of UV light for dust extinction is ^ 4 
(i.e., e ~ 0.75) for the Lyman break galaxies at z ~ 3. 
Therefore, the change of Cir as a function of z will be an 
interesting problem to which we should apply our method. 

On the UV side, GALEX data will be available in 
a few years. By applying our method to these data, we 
will be able to examine the statistical properties of UV 
extinction (or e) with the aid of the future IR data. As 
mentioned in Section 16.21 this survey will contribute to 
revealing the representative value of e for the nearby star- 
forming galaxies with the deepest UV data that we have 
ever had. We can also determine the statistical properties 
of SFR more accurately by using both UV and IR data 
(Flores et aLHTO^ Buat et al. HUU^ . 

Finally we should comment on galaxies with strong 
Lya emission, because we have assumed that all the Lya 
photons are absorbed by dust during the resonant scat- 
tering. This assumption may not be valid for objects with 
large Luv/^iR because absorption of light by dust is 
not efficient in such galaxies. Such a condition could be 
satisfied in a primeval galaxies which is little enriched 
by dust (or metals). Our assumption is only valid if 
Luv/iiR < Luy/Li^ya — 15// (this value comes from 
Starburst99 prediction). If dust extinction is inefficient, 
/ ^ 1. Therefore, if we find a galaxy with Luv/^iR ^ 15, 
it can be a candidate for a dust-deficient primeval objects 
with a conspicuous Lya emission line. 
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Appendix A: Conversion factors with a constant 
SFR over 10^ yr 

In order to see the robustness of the result against isF, 
we also construct the conversion formula for SFR for 
tsF = 10'' yr. In other words, we adopt the Starburst99 
stellar synthetic spectrum at lO'' yr with a constant star 
formation history. Since the luminosity of the Lyman con- 
tinuum photons is already stationary at 10^ yr, the quan- 
tities concerning the Lyman continuum photons and re- 
combination processes are not changed at all. We should 
however change Cir and C2ooo- If we adopt <sf — 10^ 
yr, the term "young" in should be used to indicate recent 
10^ yr instead of 10* yr as in the text. Accordingly "old" 
should be used for the age larger than 10^ yr. 



A.l. Cm 

In the main text, we have assumed the stellar synthetic 
spectrum made with Starburst99 with a constant SFR at 
the age of 10* yr. We examine the Starburst99 result for 
a constant SFR at 10^ yr with the other conditions fixed. 
As a result, we obtain ^Lyc = 0.20Lboi and Lnonion = 
O.SOLboi- Since both Li^ya and Lion are already stationary 
at 10^ yr for a constant SFH, the relation between those 
two luminosities is the same as that in the main text: 
Lhya = 0.34/LLyc- Instead of equation we obtain 

Lfl = (0.20 - 0.13/ + 0.80e) ■ (A.l) 

We also obtain the following relation instead of equation 



^™ , = 2.72 X 10-1" ^ . 
Mq yr-i Lq 



(A.2) 
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The conversion factor (eq. |SJ) is expressed as 

2 72 X 1 (^^^^(^ — n) 
^' ^) = 0.20- 0.13/ + O.8O; -(^-^^ 

The conversion factor under the dusty starburst approxi- 
mation (i.e., / = 0, e = 1, and r] = 0) becomes 

= 2.72 X 10-1" [Mq yr-i L^^] . (A.4) 
A. 2. C2000 

Because the 2000 A monochromatic luminosity continues 
to increases after 10^ yr, C2000 at 10^ yr is larger than that 
at 10^ yr. The Starburst99 result indicates that C2000 = 
3.18 X 10-^" Mq yr-i erg"! s A. 



